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Plasma  phase  transitions  (PPT)  in  dense  hydrogen  and  electron-hole  plasmas  are  investi- 
gated by  direct  path  integral  Monte  Carlo  methods  (DPIMC).  The  phase  boundary  of  the 
electron-hole  liquid  in  Germanium  is  calculated  and  is  found  to  agree  reasonably  well  with 
the  known  experimental  results.  Analogous  results  are  found  for  high-density  hydrogen.  For 
a temperature  ofT=  10, 000 A it  is  shown  that  the  internal  energy  is  lowered  due  to  droplet 
formation  for  densities  between  10 23  cm-3  and  1024  cm~3. 


1.  Path  integral  Monte 
Carlo  simulations 

All  thermodynamic  properties 
of  a two-component  plasma  are 
defined  by  the  partition  func- 
tion Z which,  for  the  case  of 
Ne  electrons  and  Np  protons, 
is  given  by  Z(Ne,Np,V,f>)  = 
Q(£;,y,  With  Q(Nc,Np,0)  = 

J dqdr  p((pr,o;P),  where 

v 


(3  = l/kpT.  The  exact  den- 

sity matrix  is,  for  a quantum 
system,  in  general,  not  known 
but  can  be  constructed  using 
a path  integral  representation 


[i,  2]  /dR<°>E.p(fl((V;/?)  = 

V 

fdRW  . . . dR^  pW  ■ pW  . ...  p<n>  x 

YJaY1p{±VKpS{<T,Pa')PPW, 
where  = p (R(l~l\  R(i^;  Afl)  = 
1)[e— whereas  A(3  = 
P/(n  + 1)  and  AA2  = 27r/i.2A/?/m„, 
a — p,e.  H is  the  Hamilton  op- 
erator, H = K + Uc,  containing 
kinetic  and  potential  energy  con- 
tributions, K and  Uc,  respectively, 
with  Uc  = tjf  + U/  + Ufr  being  the 
sum  of  the  Coulomb  potentials  be- 
tween protons  (p),  electrons  (e)  and 
electrons  and  protons  (ep).  Further, 
a comprises  all  particle  spins,  and 


the  particle  coordinates  are  denoted 
by  RW  = (gW.rW)  = (r£\r{J]), 
for  i = 1, . . .n  + 1,  R(°)  = (q,r)  = 

(Ri0),Ri0)),  and  R("+1>  = R<°> 
and  a'  — a.  This  means,  the  par- 
ticles are  represented  by  fermionic 
loops  with  the  coordinates  (beads) 
[R]  = 

where  q and  r denote  the  electron 
and  proton  coordinates,  respectively. 
The  spin  gives  rise  to  the  spin  part 
of  the  density  matrix  S,  whereas  ex- 
change effects  are  accounted  for  by 
the  permutation  operator  P,  which 
acts  on  the  electron  coordinates  and 
spin  projections,  and  the  sum  over 
the  permutations  with  parity  Kp. 
To  compute  thermodynamic  func- 
tions, the  logarithm  of  the  partition 
function  has  to  be  differentiated 
with  respect  to  thermodynamic 
variables,  so  for  internal  energy  E 
we  have  (3E  = —fidlnQ/dfi 
2.  Numerical  Results 
Since  the  PPT  in  dense  hydrogen 
is  still  hypothetical  and  has  not 
been  observed  experimentally,  it  is 
reasonable  to  look  for  other  sys- 
tems where  similar  conditions  exist. 
A suitable  example  is  electron-hole 
plasma  in  low-temperature  semicon- 
ductors, for  which  droplet  formation 
is  well  established  and  observed  ex- 


3 


perimentally  three  decades  ago  [3]. 
We,  therefore,  performed  DPIMC 
simulations  for  electron  hole  plas- 
mas. Below  the  critical  tempera- 
ture the  simulations  exhibit  anoma- 
lously large  fluctuations  and  an  un- 
stable behavior  of  the  pressure.  The 
e-h-plasma  is  found  to  phase  sepa- 
rate and  form  large  droplets.  The 
phase  boundary  of  the  electron-hole 
liquid  (e-h-droplets)  in  Germanium 
obtained  by  our  DPIMC  method  is 
presented  in  Fig.  1 together  with  the 
experimental  data.  We  observe  good 
agreement.  Deviations  may  be  con- 
nected with  complex  band  structure 
of  Germanium  approximated  in  our 
simulations  by  a two-band  parabolic 
mass  model. 


Figure  1:  Phase  boundary  of  the 
electron-hole  liquid  in  bulk  Germa- 
nium. Experiment  - [3].  Tempera- 
ture is  presented  in  units  of  the  ex- 
citon  binding  energy. 

Fig.2  presents  results  of  our  calcu- 
lations for  pressure  and  energy  of 
hydrogen  at  T = 10, 000 R,  which 
is  well  below  the  critical  point  of 
the  PPT  predicted  by  chemical  mod- 
els (around  T = 15,000K).  In  our 
calculations  between  1022cm~3  and 
1024cm-3  calculated  pressure  be- 
comes negative.  We  find  also  that 
in  this  region  plasma  energy  is  sys- 
tematically lower  than  the  RPIMC 
results  [4],  If  average  distance  be- 


Figure 2:  Pressure  and  energy  of  hy- 
drogen for  T = 10,000/s:,  DFT  - [5], 
RPIMC  - [4]. 


tween  plasma  particles  is  of  the  or- 
der of  the  size  of  a hydrogen  molecule 
the  homogeneous  plasma  state  be- 
comes unstable,  and  many-particle 
clusters  appear.  Results  of  one  in- 
dependent well  tested  method  based 
on  density  functional  theory  (DFT) 
are  presented  on  Fig.  2,  where  PPT 
was  obtained  at  smaller  density  [5]. 
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